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Abstract 

We are looking for the agent-based treatment of the financial markets considering necessity to build bridges between 
microscopic, agent based, and macroscopic, phenomenological modeling. The acknowledgment that agent-based 
modeling framework, which may provide qualitative and quantitative understanding of the financial markets, is very 
ambiguous emphasizes the exceptional value of well defined analytically tractable agent systems. Herding as one of the 
behavior peculiarities considered in the behavioral finance is the main property of the agent interactions we deal with in 
this contribution. Looking for the consentaneous agent-based and macroscopic approach we combine two origins of the 
noise: exogenous one, related to the information flow, and endogenous one, arising form the complex stochastic dynamics 
of agents. As a result we propose a three state agent-based herding model of the financial markets. From this agent-based 
model we derive a set of stochastic differential equations, which describes underlying macroscopic dynamics of agent 
population and log price in the financial markets. The obtained solution is then subjected to the exogenous noise, which 
shapes instantaneous return fluctuations. We test both Gaussian and q-Gaussian noise as a source of the short term 
fluctuations. The resulting model of the return in the financial markets with the same set of parameters reproduces 
empirical probability and spectral densities of absolute return observed in New York, Warsaw and NASDAQ OMX Vilnius 
Stock Exchanges. Our result confirms the prevalent idea in behavioral finance that herding interactions may be dominant 
over agent rationality and contribute towards bubble formation. 
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Introduction 

Statistical physics has got the edge over socio-economic sciences 
in the understanding of complex systems [1-7]. This happened 
due to the fact that physicists were able to start from the 
understanding of simple phenomena via simple models and later 
built the complexity up together with the increasing complexity of 
the considered phenomena. On the other hand socio-economic 
sciences had to face complexity right from the start as socio- 
economic systems are in no way simple systems - they are 
intrinsically complex at many different levels at the same time. 
Financial markets are one of the most interesting examples of such 
complex systems. Unlike in physics we have no direct way to gain 
insights into the nature of microscopic interactions in financial 
markets, thus the understanding of the financial market fluctua- 
tions may become rather limited and very ambiguous. Yet the 
understanding might be improved indirectly through the further 
development of the complex systems approach [8]. First of all, 
currently there are huge amounts of the available empirical data, 
which itself is attracting representatives of the experimental 
sciences [9]. Also there is an agent-based modeling framework, 
which may provide qualitative and quantitative understanding of 
the financial markets. The intense applications of these ideas is still 
ongoing [9-1 1] and the challenge is still open. 



Agent-based modeling has become one of the key tools, which 
could improve the understanding of the financial markets as well 
as lead to the potential applications [12-17]. Currently there are 
many differing agent-based approaches in the modeling of the 
financial markets. Some of them aim to be as realistic as possible, 
yet they usually end up being too complex to posses analytical 
treatment. One of the most prominent examples of these kind of 
models is the Lux-Marchesi model [18]. The more recent 
approaches in the similar direction consider modeling order books 
[19-21]. Other approaches, on the other hand, aim to capture the 
most general properties of the many complex socio-economic 
systems (some of the examples include [22,23]). Though there are 
also some interesting approaches which combine realism and 
analytical tractability, e.g. Feng el al. [24] have used both 
empirical data and trader survey data to construct agent-based 
and stochastic model for the financial market. Looking for the 
ideal agent-based approach we would consider as a primary 
necessity to build bridges between microscopic, agent based, and 
macroscopic, phenomenological, modeling [25]. Following this 
trace of thought it would be rational to combine two origins of the 
noise: exogenous one, related to the information flow, and 
endogenous one, arising form the complex stochastic dynamics 
of agents. Such integral view of the financial markets can be 
achieved only with very simple zero-intelligence agent-based 
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models and macroscopic, phenomenological, approaches incor- 
porating external information flow. This is the main idea of our 
present consideration of the financial markets. 

The expected properties of such model lead us to the return 
fluctuations, characterized by the power law distributions and the 
power law autocorrelations of absolute return considered in [26- 
32]. We investigate an agent-based herding model of the financial 
markets, which proves to be rather realistic and also simple 
enough to be analytically tractable [33,34] . Namely we consider a 
three agent states' model [34] and incorporate it into the standard 
model of the stock price described by the geometric Brownian 
motion or into process with statistical feedback [35], exhibiting 
Tsallis statistics. We find that the improved three state agent-based 
herding model reproduces the power law statistics observed in the 
empirical data extracted from the NYSE Trades and Quotes 
database, Warsaw Stock Exchange and NASDAQ, OMX Vilnius 
Stock Exchange. 

We start by discussing the possible alternatives in macroscopic 
and phenomenological modeling providing some insight into the 
possible connection to the agent-based microscopic approach. 
Next we develop the microscopic approach by defining the 
herding interactions between three agent groups and incorporate it 
into a consistent model of the financial markets. Further we couple 
the endogenous fluctuations of the agent system with the 
exogenous information flow noise incorporated in macroscopic 
approach and provide detailed comparison with the empirical 
data. Finally we discuss the obtained results in the context of the 
proposed double stochastic model of the return in the financial 
markets. 

Methods 

Macroscopic and phenomenological versus microscopic 
and agent-based treatment of the financial markets 

It is the natural peculiarity of the social systems to be treated 
first of all from the macroscopic and phenomenological point of 
view. In contrast to the natural sciences microscopic treatment of 
the social systems is ambiguous and hardly can be considered as a 
starting point for the consistent modeling. The complexity of 
human behavior leaves us without any opportunity to consider 
human agent in action as a determined dynamic trajectory. The 
financial markets as an example of the social behavior first of all 
are considered as a macroscopic system exhibiting stochastic 
movement of the variables such as asset price, trading volume or 
return [10]. Despite lack of knowledge regarding microscopic 
background of the financial systems there is considerable progress 
in stochastic modeling producing very practical applications 
[36,37]. The standard model of stock prices, S(t), referred to as 
geometric Brownian process, is widely accepted in financial 
analysis 

dS = n(t)Sdt + a(t)SdW. (1) 



processes [38-40] as well as for the stochastic modeling of volatility 
u(t) [37]. 

We acknowledge this phenomenological approach as a good 
starting point for the macroscopic financial market description 
incorporating external information flow noise W and we will go 
further by modeling volatility a(t) as an outcome of some agent- 
based herding model. The main purpose of this approach is to 
demonstrate how sophisticated statistical features of the financial 
markets can be reproduced by combining endogenous and 
exogenous stochasticity. 

Here we consider only the most simple case, when a(t) 
fluctuations are slow in comparison with external noise W . In 

S(t+ T) 

such case the return, r t (T) = In — — — — , in the time period T can 

S(t) 

be written as a solution of Eq. (1) 

r t {T)= (pit)- \<j(t) 2 ^j T + a(t)W(T). (2) 

This equation defines instantaneous return fluctuations as a 

1 , 

Gaussian random variable with mean (fi— - a )T and variance 

a 2 T. Let us exclude here from the consideration long term price 

1 , 

movements defined by the mean (p. — — ff ) T as we will define the 

dynamics of price from microscopic agent-based part of model. 
This assumption means that we take from phenomenological 
model only the general idea how to combine exogenous and 
endogenous noise. Then Eq. (2) simplifies to the instantaneous 
Gaussian fluctuations N[0,a(t) 2 T] with zero mean and variance 
o(t?T. 

In [41,42], while relying on the empirical analysis, we have 
assumed that the return, r,(T), fluctuates as instantaneous q- 
Gaussian noise N q [ro(x),X\ with some power-law exponent 
2 

A = = 5, and driven by some stochastic process x(t) defining 

second parameter of fluctuations ro(x). ro(x) was introduced as a 
linear function of absolute return moving average \x\ calculated 
from some nonlinear stochastic model [42] 

r 0 (x) = b + a\x\, (3) 

where parameter b serves as a time scale of exogenous noise and - 

quantifies the relative input of exogenous noise in comparison with 
endogenous one described by \x(t)\. 

A more solid background for this kind of approach can be found 
in the work by L. Borland [43]. The idea to replace geometric 
Brownian process of market price by process with statistical 
feedback [35] leads to the equation of return r,(T) as function of 
time interval T given by 



In the above Wiener process W can be considered as an 
external information flow noise while a(t) accounts for the 
stochastic volatility. Though one must consider the model of 
stock prices following geometric Brownian motion as a hypothesis 
which has to be checked critically, this serves as a background for 
many empirical studies and further econometric financial market 
model developments. Acknowledgment that analysis taking pi(t) 
and o( t) constant have a finite horizon of application has become 
an important motivation for the study of the ARCH and GARCH 



&r t {T) = o(t)&Q.T, (4) 

where Qj evolves according to the statistical feedback process [35] 

dn r =p(n,rr?W r , (5) 

and P(Q,T) satisfies the nonlinear Fokker-Planck equation 
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-^-P(Q,T) = -^rP 2 -i(Q.,T). 

St y ' an 2 



(6) 



The explicit solution for P in the region of q values 1 <q< 5/3 
can be written as one of the Tsallis distributions [44] 



p g (n,T)= ^{i-y(r)(i- 9 )n 2 }i\ 

where C(T) and y(T) are as follows 
3-q 



(7) 



C(T) = 



VST 



2(9-1) 



i 



2 

3-? 



[(2-q)Ci-q)TY-9, 



y(T) = 



q-l 



3-g 
2(9-1) 



2(?-l) 
3-9 



g-1 

(9-1)3-9 



[(2-9X3-4)713-? 



(8) 



Assuming ct(?) as slow stochastic process in comparison with Qt 
from Eq. (4) one gets that r,(T) = a(t)Q T . This sets PDF for r,(T) 
the same as for Qt Eq. (7), one just has to replace Q and, C(T) 



defined in Eqs. (8) by 



r,(T) 



and a(t)C(T) accordingly. This gives a 



Tsallis distribution for r,(T) as 



P q (r,{T))- 



1 



a(t)C(T) 



!-(!-<?) 



r,(Tf 



0-q)a(t) 2 a q {T) 2 



l 

1=9 



(9) 



where a q {T) as new parameter related to previous one f(T) can 
be written as 



o q (T) = 



3-g 
3-«f V2(g-1). 



(9-1) 



q- 



y/n(X-l)- 



q-l 



[(2- ? )r]3-9= (io) 



i-i 



U-2), 



(11) 



asset deviates from the fundamental value. These deviations can 



be accounted as p(t) = In 



Pit) 

Pf(ty 



Further in this contribution we 



will assume that volatility a(t) is defined by \p(t)\ through the 
linear relation 



<r(t) = b{\+a\p(t)\), 



(12) 



where parameter b serves as a scale of exogenous noise and a 
quantifies the relative input of endogenous noise. Both parameters 
a and b have to be defined from empirical data. To complete the 
model we have to propose agent-based consideration of log price 
p(t). In the following section we present the three state herding 
model giving stochastic equations for the log price p(t). 

The three state herding model as a source for the 
endogenous stochastic dynamics 

Having discussed a macroscopic view of the financial fluctua- 
tions we now switch to the microscopic consideration of the 
endogenous fluctuations. Let us derive the system of stochastic 
differential equations defining the endogenous log-price fluctua- 
tions from a setup of appropriate agent groups composition. We 
consider a system of N heterogenous agents - market traders 
continually changing their trading strategies between three 
possible choices: fundamentalists, chartists optimists and chartists 
pessimists. We further develop this commonly used agent group 
setup [24] by considering all transitions between three agent states 
to be a result of binary herding interactions between agents during 
their market transactions. 

Fundamentalists are traders with fundamental understanding of 
the true stock value, which is commonly quantified as the stock's 
fundamental price, Pf(t). We exclude from our consideration any 
movements of the fundamental price. The assumption of constant 
fundamental price means that we further will consider price 
fluctuations around its fundamental value. Taking into account a 
long-term rational expectations of the fundamentalists their excess 
demand, EDf(t), might be assumed to be given by [45] 



ED f (t) = N f (t)\n 



pay 



(13) 



where N/(t) is a number of the fundamentalists inside the market 
and P(t) is a current market price. The rationality of fundamental 
traders keeps asset price around its fundamental value as they sell 
when Pf<P{t), and buy when Pf>P(t). 

Short term speculations cause unpredictable price movements. 
There is a huge variation of speculative trading strategies, so it is 
rational, from the statistical physics point of view, to consider these 
variations as statistically irrelevant. We make only two distinctions 
between chartists: optimists suggest to buy and pessimist suggest to 
sell at a given moment. Thus the excess demand of the chartist 
traders, ED c (t), can be written as: 



Now we are prepared to combine two phenomenological 
approaches introduced by Eqs. (2) and (4) with agent-based 
endogenous three state herding model. o(i) serves as a measure of 
system volatility in both of the phenomenological approaches. It is 
reasonable to assume that financial market is in the lowest level of 
possible volatility when assets market value P(t) is equal to the it's 
fundamental value Pf(t), lets define it as constant a(t) = b. 
Volatility of financial system increases when market value of the 



ED c (t) = r[N 0 (t)-N p (t)], 



(14) 



where r is a relative impact factor of the chartist trader and further 
will be integrated into a certain empirical parameter, N 0 and N p 
are the total numbers of optimists and pessimists respectively. So, 
as you can see, we replace a big variety of "rational" chartist 
trading strategies by herding kinetics between just two options: buy 
or sell. 
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The proposed system of heterogenous agents defines asset price 
movement by applying the Walrasian scenario. A fair price reflects 
the current supply and demand and the Walrasian scenario in its 
contemporary form may be written as 



variable rif changes slowly in comparison with n p or n 0 is helpful 
here. This enables to consider n p dynamics as one dimensional 
process as well. Let us write the transition rates for n p in the same 
way as in Eq. (18) 



1 dp(Q 
pN dt 



= -n f (t)p(t) + r[n 0 (t)-n p (t)], 



(15) 



here p is a speed of the price adjustment, N a total number of 

traders in the market, p(t) = In — ^ and = — By assuming 

Pf N 

that the number of traders in the market is large, JV— >oo, one 

obtains: 



P(t)- 



_n 0 (t)-n p (t) 
n f (t) 



(16) 



n po (n p ) = n p ^ + Hh( 1 - n f - n p )^ , 

x op (n p ) = (l-n f -n p )^-+Hhn p y (20) 



As in other similar cases [33] these one step transitions lead to 
the SDE for n„ 



dn p = (1 — iif — 2n p )a cc dt + ^j2Hhn p {\—nf — n p )dW p . (21) 



Stochastic dynamics of the proposed agent-based system is 
defined by the occupations of the three agent states: 



N ' 



N, 



p 
N ' 



No 
N ' 



(17) 



One can model the evolution of occupations as a Markov chain 
with some reasonable assumptions for the sake of simplicity. There 
are six one agent transition possibilities in three group setup, see 
Fig. 1 and [34]. Few assumptions are natural for the financial 
markets as there is some symmetry. With the notation of agent 
transitions from state i to j as subscripts to any parameter Ay, 
where i and j take values from the set {f,c,o,p} = 
{fundamentalists, all chartists, optimists, pessimists}, we will 
use following assumptions, for the rates of spontaneous transitions: 
Oop = Opo = o C c, Of P = Of 0 = af c /2, <T p f = a 0 f = a c f, and for the 
herding transitions: hf p = h p f = h/ 0 = h 0 f = h. Finally it is reasonable 
to assume that transitions between chartist states are much faster 
than between chartists and fundamentalists h op = h po = Hh, H»\, 
a cc » a c f, o C c » o/c- Taking into account the restraint Nf + N C = N 
and having in mind that transitions f— >o axe equivalent to f-*p, 
one can can write one step herding transition rates between / and 
c groups for given rif = Nf/N and n c = N c /N as [33] 



n fc (n f ) = n f ^ +h(l- n f ) 
n c f(n f ) = ( 1 - nf) + to/] . 



(18) 



Here transition rates have the same form as in Kirman's 
herding model [46]. Using Eq. (18) one can write the Master 
equation for PDF P(tif,i) and derive Fokker- Planck equation in 
the limit JV— >go [33]. One dimensional Fokker-Planck equations 
has its equivalent stochastic differential equation which for the nf 
can be written as 



dn f = [( 1 - n f )o c f - n f o fc ] dt + ^2hn f {\ -n f )dW f . (19) 



The next step in this approach is to define dynamics of n p under 
restraint nf + n p + n 0 = \. An adiabatic approximation assuming 



Equations (19) and (21) form a system of coupled SDEs and 
define the agent population dynamics in three state agent model 
taking into account the previous assumptions. It is possible to 
rewrite Eq. (21) in the form without direct dependance on My by 

introducing another variable £(t) = \ average mood of 

8 n 0 (t)+n p (t)' 8 

the chartists, instead of n p . Such variable substitution makes 

second SDE independent of the first one 



d£ = - 2o c Mt + \j2Hh(\-ii 2 )d W t . 



(22) 



Equations (19) and (22) are independent and define occupation 
dynamics of the same three state agent-based herding model. 

In the previous work we considered generalization of the 
herding model introducing variable interevent time T, see [33,34]. 
It is a natural feedback of a macroscopic state on the microscopic 
behavior, activity of agents. Such feedback is an empirically 
defined phenomena and can be quantified through the relation of 
trading volume with return [47]. We introduce this feedback into 
proposed three state herding model as a trading activity, rate of 
1 

transactions, — — defined as 



l+a 



\—Hf 



= (l4-4p(0ir, (23) 



where empirical parameter a is the same as in Eq. (12). Such power- 
law behavior is consistent with empirical data quantifying relation of 
short term return r with trading volume V, V(r)~r a , where a~2, 
[48-50]. Notice that for the geometric Brownian motion model of 
stock price the return as increment of log-price is proportional to the 
log-price. This is the another reason together with search of 
simplicity we use log-price instead of return in Eq. (23). 

Taking into account the discussed feedback mechanism, intro- 
ducing scaled time, t s = ht, and appropriately redefining model 
parameters: e c f = a c f/h, £f c = Of c /h, e cc = a cc /(Hh), we are able to 
rewrite stochastic differential equations of endogenous model as 



dnf 



(l-n f )e cf 



dt s + 



2nf{\—nf) 
-c(n f ,0 



dW s 



(24) 
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Figure 1. Schematic representation of the three state herding model, where relevant parameters are shown. The arrows point in the 
directions of the possible transitions, each of the transitions pairs is modeled using original Kirman's model. 
doi:1 0.1 371 /journal.pone.01 02201 .g001 



t(«/,<D ' 



dL 



277(1 



-<*) 



(25) 



Ci +i =£r 



2hH£ cc £i 



At t - 



2/i77(l 



(28) 



Then log-price defined in Eq. (16) can be expressed through the 
stochastic variables of the model iif and £ 



7*0=^(0. 

«/(0 



(26) 



In Eq. (26) and further we omit parameter f as consider it 
integrated into parameter a. 

This concludes the definition of consentaneous agent-based and 
stochastic model of the financial markets as Eq. (26) defines joint 
endogenous and exogenous volatility a(t) introduced by Eq. (12). 

Numerical algorithm 

We solve Eqs. (24) and (25) by using Euler-Maruyama method 
[51] with variable time step, A?,. Namely we have transformed a 
set of stochastic differential equations into a set of difference 
equations: 



= Xj + h 



-XiSfi 



At,- 



2hxi{ \ — Xj) 



i(xi,£i) 



A^ u , (27) 



t i+ i = ti+Ati, At r - 



K 2 T(Xi,ii) 

h{ 1 + e cf +£f c + 77(1 + 2e cc )) ' 



(29) 



In the above we have changed the notation, Xt = rifj, to improve 
readability of the difference equations. While Cjj stands for the 
uncorrelated normally distributed random variables with zero 
mean and unit variance. Also note that in the difference equations 
above we have introduced additional parameter K, which is 
responsible for the precision of the numerical results. The smaller 
K value gets, the more precise numerical simulations are, but the 
longer computation time grows. During the numerical simulations 
we found K = 0.03 to be the optimal value precision-wise and time- 
wise. 

In order to keep «/ and £ well defined we introduce absorbing 
boundaries near the edges of the intervals in which these variables 
are defined. Namely, we require that each iif belongs to [5,1 —6] 
and £ belongs to [— 1 + <5,1 — 8}. In other words before using 
and obtained from the Eqs. (27) and (28), we apply min( . 
and max( . . . ) functions on them: 



= min(max(jt,- + 1 ,3), 1 — 5), 



(30) 
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<f,-+i=min(max({ i+1 ,-l+<5),l-(5). (31) 

The new i and £',+1 are certainly well defined and may be 
used in further simulations. In our simulations we have used 

£>= 10 6 . 

Next step in our numerical simulation is to obtain time series of 
Pi, discretized at one minute time periods, and apply ^-Gaussian, 

r(t,T) = b(l+a\ 1 ^i t \y q (T)N q (0,l), (32) 
or, in certain cases, Gaussian, 

r(t,T) = b(^+a\ ] ^^VfN(0,l), (33) 

noise on it. Afterwards we normalize the obtained one minute 
return time series so in the end it would have unit variance. 

Finally we add as many one minute returns as we need to obtain 
the final return time series at a desired discretization intervals. In 
the second section we analyze empirical data from New York, 
Warsaw and NASDAQ OMX Vilnius Stock Exchanges and 
reproduce it's statistical properties in the numerical simulation. 

Results and Discussion 

Comparison of model simulation and empirical data 

Now we will adjust model parameters to reproduce empirical 
data of the return in three different markets. The model return 
r t {T) in the time interval T can be written as 

r(t,T) = a(t)VfN(0,l), (34) 
for the Gaussian external noise, see Eq. (2), and 

r{t,T) = N q {Q,o{i)o q {T)), (35) 

for the q-Gaussian one, see Eq.(9). Here N(Q,\) denotes normally 
distributed random variable with zero mean and unit variance, 
and N q (0,a(t)a q (T)) denotes Tsallis random variable distributed 
as defined by Eq.(9). We choose T=\ minute as primary 
sufficiently short time interval where ff(t) fluctuations are 
negligible and solve equations (24) and (25) numerically in 
successive time intervals to get 1 minute time series of p(t), Eq. 
(26). From the price time series one can produce time series for the 
volatilities a(t) and returns r{t,T) using Eq. (34) or Eq. (35). 

We compare model return series with empirical return time 
series extracted from high frequency trading data on New York, 
Warsaw, and NASDAQ, OMX Vilnius Stock Exchanges. These 
series were transformed into successive sequences of empirical 1 
minute returns. Produced empirical return series were normalized 
by standard deviation calculated on the entire time sequence of 
selected stock. For this comparison with empirical data we select 
only a few stocks from each stock exchange, which have more or 
less constant long term average trading activity to avoid 
considerable input of possible trend into time series statistics. 
From NYSE data we have selected stocks: BMY, GM, MO, T, 
traded for 27 months from January, 2005. From Warsaw SE 
stocks: TPSA, KGHM, traded from Novemver 2000 to January 
2014, and PZU traded from May 2010 to January 2014. From 



NASDAQ OMX Vilnius data stocks: APG1L, IVL1L, PTR1L, 
SAB1L, TEOIL, traded from May 2005 to December 2013. We 
do not extend this analysis into more wide representation of stocks 
as only few stocks from NASDAQ, OMX Vilnius are liquid 
enough for such analysis. Comparable demonstration of general 
features across markets is the main purpose of this consideration. 
In every stock group series of different stocks are considered as 
separate realizations of the same stochastic process and so we 
present empirical statistical information as average over realiza- 
tions (stocks). 

In order to compare scaling of the statistics with increasing time 
interval T of return definition we just sum 1 minute returns of 
model and empirical time series in successive intervals of 3 
minutes, of 10 minutes or 30 minutes. For each stock exchange 
considered and each time interval T we calculate absolute return 
probability density function (PDF) and power spectral density 
(PSD). Both PDF and PSD are obtained by averaging over stocks 
in the considered stock exchanges. Results are presented in four 
figures: Fig. 2 - for NYSE data; Fig. 3 - for NYSE data with 
Gaussian noise; Fig. 4 - for Warsaw Stock Exchange data; Fig. 5 - 
for NASDAQ OMX Vilnius data. In figures: Fig. 2, Fig. 4, Fig. 5 
empirical data (red lines) are compared with model statistics (black 
line) calculated with the same choice of parameters: £ c f = 0.l, 
e fc = 3, e cc = 3, 77 = 300, A=10- 8 j- 1 , a = 0.5, b = \, A =4, <x = 2. 

From our point of view the achieved coincidence of empirical 
and model statistics is better than expected. PDF and PSD 
coincide almost for all markets and all time T scales. q-Gaussian 
noise suites much better for this model, compare Fig. 2 calculated 
with q-Gaussian noise and Fig. 3 calculated with Gaussian one. 
The use of q-Gaussian noise can be confirmed by more detailed 
study of empirical data as well. 

There are some observed discrepancies of empirical and model 
results, which have reasonable explanations. For example, spikes 
observed in empirical PSD are related with one trading day 
seasonality, this is not included in the presented model and not 
observed in the model PSD. Thanks to reviewer, who has 
stimulated us to think how could be empirically observed intra-day 
pattern of return volatility and trading activity introduced into our 
model. The answer is not so simple as one could expect because 
we use simplified, assumed as only statistically valid, relation 
between log price p(i) and trading activity, see Eq.(23). The 
proposed model does not include trading activity as independent 
dynamic variable and this makes introduction of discussed 
seasonality regarding trading activity not straightforward. Never- 
theless, this intraday seasonality first of all has to be reflected in 
volatility pattern of total return fluctuations quantified by 
parameter b. Let us to replace constant value of b = 1 by some 
intra-day exponential variation 

b(t)= exp [-(r— 195) 2 /w 2 ] +0.5, (36) 

where w = 20 quantifies in minutes the width of b variation and 
time t is closed in the circle of total duration of NYSE trading day 
equal to 390 minutes, note ( = 0 or ^ = 390 corresponds to the 
middle point of trading day. In Fig. 6 we present results of such 
numerical experiment, where power spectral density of absolute 
return for NYSE stocks is compared with model accounting for b 
seasonality, Eq. (36). As one can see, such a simple introduction of 
volatility seasonality reproduces empirical resonance structure of 
absolute return very well. 

For such emerging markets as Baltic region the intensity of stock 
trading is up to 50 time lower than in NYSE. This makes a 
considerable impact on return in such short time interval as 1 
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minute PDF as for over 90% of time intervals no transactions are 
registered. Consequently, one observes considerable discrepancy 
in PDF for low return values as main part of returns are equal to 
zero. Nevertheless, it is worth to notice that power law part of PDF 



and PSD is in pretty good agreement with model one. Mentioned 
discrepancy disappears with increasing T or for the markets with a 
more intensive trading. Some PDF discrepancy for low returns 
values is observed for NYSE and W arsaw stocks as well and this is 
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Figure 3. Probability density function (PDF) and power spectral density (PSD) of return with Gaussian noise for NYSE stocks: BMY, 
GWI, MO, T. Empirical (red) and model (black) PDF (first column) and PSD (second column), (a) and (b) - 1 minute; (c) and (d) - 3 minutes; (e) and (f) - 
10 minutes; (g) and (h) - 30 minutes. Model parameters are as follows: e c j =0.1, £/ c = 3, e cr = 3, H = 300, /) = 10~ s .s~', a = 0.5, 6 = 1, a=2, 
doi:1 0.1 371 /journal.pone.01 02201 .g003 



related with some impact of discrete nature of price values 
measured in cents, namely on the discrete tick size. Certainly, 
prices are smooth in the model. 



Finally we can conclude that three state herding model of 
absolute return in financial markets works very well and explains 
general origins of power law statistics for very different markets 
starting from most developed to emerging. 
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Figure 4. Probability density function (PDF) and power spectral density (PSD) of return with (/-Gaussian noise for Warsaw Stock 
Exchange stocks: KGHM, PZU, TPSA. Empirical (red) and model (black) PDF (firs column) and PSD (second column), (a) and (b) - 1 minute; (c) and 
(d) - 3 minutes; (e) and (f) - 10 minutes; (g) and (h) - 30 minutes. Model parameters are the same as in Fig. 2. 
doi:1 0.1 371 /journal.pone.01 02201 .g004 



Conclusions 

The proposed consentaneous agent-based and stochastic model 
of the financial markets is a result of our previous research in 



stochastic modeling, see references in [25,41], and agent-based 
modeling of herding interaction [33,34]. The idea to combine two 
approaches comes from the properties of nonlinear stochastic 
differential equations generating power law statistics [41,52] and 
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Figure 5. Probability density function (PDF) and power spectral density (PSD) of return with (/-Gaussian noise for NASDAQ OMX 
Vilnius stocks: APG1 L, IVL1 L, PTR1 L, SAB1 L, TEOI L. Empirical (red) and model (black) PDF (first column) and PSD (second column), (a) and (b) - 
1 minute; (c) and (d) - 3 minutes; (e) and (f) - 10 minutes; (g) and (h) - 30 minutes. Model parameters are the same as in Fig. 2. 
doi:1 0.1 371 /journal.pone.01 02201 .g005 



possibility to derive these equations as macroscopic outcome of different exponents distinguishes this approach as compromise of 

microscopic herding model [33]. The possibility to reproduce sophistication and simplicity. Nevertheless, the main value of this 

power law statistical features of absolute return observed in the model comes form its very clear microscopic background of 

financial markets in such details as power spectrum with two herding interactions between agents. Very general idea, which has 
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Figure 6. Power spectral density (PSD) of return with (/-Gaussian noise and seasonality accounted for NYSE stocks: BMY, GM, MO, T. 

Empirical (red) and model (black) PSD accounting for parameter's b seasonality Eq. (36). (a) - 1 minute; (b) - 3 minutes; (c) - 10 minutes; (d) - 30 
minutes. Other model parameters are the same as in Fig. 2. 
doi:1 0.1 371 /journal.pone.01 02201 .g006 



roots in the entomological studies of ant colonies [46], can be 
adopted to build system of three agent groups acting in the 
financial markets. Further complexity of the behavioral aspects of 
agents can be treated as irrelevant for the macroscopic outcome of 
this complex system as herding alon reproduces statistical 
properties. 

It is worth to notice that some feedback of observed return 
increases trading activity in the market and the degree of 
nonlinearity, which we account by empirical value of the exponent 
a = 2, [48-50] . Other model parameters are less grounded by the 
empirical market analysis and are defined by adjusting proposed 
model to the here presented statistical properties of the considered 
financial markets. Parameters e c t = 0.1 and et c = 3 define the most 
fundamental tradeoff between fundamental and speculative 
behavior of agents in financial markets. It is clear from the 
defined values of the parameters, exhibiting strong asymmetry 
between fundamentalism and chartism, that bubbles in the 
financial markets can be explained just by the disappearance of 
the fundamental trading behavior and markets become stable 
when considerable part of agents return to the trading according 
to the fundamental values. The tradeoff between fundamental and 
speculative behavior is the slowest process in the model defined by 
adjusted parameter h = 10~ 8 .y~'. It looks like that this long term 
herding process has global nature and impacts all markets in the 
similar way. 

Much more rapid process defined by parameter H = 300 
describes tradeoff between choices to sell or buy and probably is 
more related with local stock price dynamics. Nevertheless, the 
same value of H appropriate for all markets and all stocks makes 
some surprise and probably is related with universal behavior of 
power spectral density. Notice that value e cc = 3 is higher than 1 



making these rapid fluctuations symmetric and localized around 
equilibrium £ = 0. 

From our point of view this study based on the concepts of 
statistical physics contributes to the behavior finance supporting 
the general idea of market inefficiency [53,54] as exhibits possible 
dominance of herding interactions over agent's rationality. In 
other words, rationality is too heterogeneous to resist herding 
tendencies. This rises the question whether markets are able to 
determine true values of assets or one needs other more 
fundamental sources of economic information. Fortunately the 
tendencies of imitation open the new possibility that herding itself 
can be used to stabilize unwanted fluctuations of the financial 
markets [55]. 

Supporting Information 

Dataset SI Compressed set of time series for one 
minute empirical returns. Time series for each stock are 
given as separate, txt files, file names are self-explanatory. 
(ZIP) 
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